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ABSTRACT 


The unique properties of the temperature and salinity profiles for polar 
oceans are critical for high-latitude mixed layer thermodynamics. In the Polar regions the 
water column is coldest and freshest at the surface where ice may be present. This 
density structure often leads to entrainment and affects both the mixed layer depth and 
the ice thickness. 

Thermobaricity, the combined dependence of seawater thermal expansion 
on temperature and pressure, magnifies the buoyancy flux associated with mixed layer 
convection. When thermobaricity amplifies entrainment so that the heat into the mixed 
layer is greater than the heat leaving the water column, the mixed layer warms and any 
existing ice begins to melt. Similarly, if the heat entrained is less than the heat leaving 
the column, the mixed layer cools and freezing occurs at the surface. In the former 
situation a polynya, or region of no ice surrounded by ice coverage, may form. 

A one-dimensional vertical model is built, and trial cases are run to show 
the intricate relationships that govern the heat and salt fluxes and subsequent ice 
thickness. The model shows the importance of thermobaricity to the air-sea-ice 
interactions. It also offers significant insight into how relatively constant atmospheric 
forcing can lead to polynya-like conditions. 
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I. INTRODUCTION 


Freezing and melting of iee in the Aretie and Antaretie environments are 
determined by the thermodynamie fluxes among air, sea-iee, and oeeanie mixed layer. 
The energy exehanged among these three domains is in a delieate balanee that may 
theoretieally lead to an iee-free Aretie basin or another iee age in the limits. This paper 
will foeus primarily on high-latitude iee-mixed layer dynamies in both the northern and 
southern hemispheres. Although the vertieal energy fluxes within these regions are 
governed by similar underlying physics one dimensionally, the two polar regions have 
differences. One phenomenon possible in all polar regions is the polynya. 

Polynyas may form as a local extreme state on a small scale, without major large- 
scale or basin changes. As local extremes, polynyas are a potential laboratory for climate 
change. Thus, part of the of the impetus for this research is stated clearly by Smith et al. 
(1990), “At predictable, recurrent locations throughout the polar regions, there are 
oceanic areas which remain either partially or totally ice free at times and under 
climatological conditions where we would expect the water to be ice covered.” 
Conversely, ice may also form in these regions at times when climatological conditions 
might warrant no-ice conditions. Simply put, there is more to ice formation and melting 
than climatological conditions; the understanding of the dynamics and thermodynamics 
of the water column below and coupled with the ice needs to be explicitly included into 
prediction models. 

This work has many implications for the United States Navy both for now and in 
the future. For the present time a better understanding of the physics of the mixed layer 
will improve ice thickness and ice extent forecast models like PIPS (Preller et al, 2002). 
The output from such a model could be vital for submarine and or shipping operations in 
high-latitude regions. Furthermore, information gained from these models could improve 
mission effectiveness and efficiency while preventing the wasting of valuable resources 
like equipment, time, and money. For example, if a scenario were developed where a P3- 
Orion was required to lay sonobuoys in the Barents Sea it would be very useful in the 
planning stages of this mission to have information regarding the ice coverage. If the 
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region is ice covered, an ice penetrating sonobuoy would be very useful assuming the 
thickness of the ice does not exceed the burying capabilities of the sonobuoy. However, 
if the region were either ice-free or the ice were too thick then the employment of this 
type of sonobuoy would be a waste of time and equipment. Similarly, the use of a 
conventional sonobuoy in an ice-covered region would be useless. Either way real-time 
knowledge about the ice thickness and extent would be invaluable for the planning of this 
mission. 

The long-term importance of this study, for the Navy, may go beyond forecast 
capabilities and mission effectiveness to include macro-scale circulation and deep water 
formation. The same physical processes that are explored in this work as mechanisms for 
ice formation and melting also seemingly play a significant role in taking oxygen and 
nutrients from the surface and reintroducing them back into the water column. When 
these nutrients are returned deep in the water column they flow with the denser water as it 
spreads around the globe replenishing the world’s oceans. Understanding completely 
how this process works and the resulting ramifications for oceanic and atmospheric 
forecasts could provide a decisive tactical advantage. 

Chapter II explores the regional descriptive oceanography of the high latitudes. 
Climatology, basin circulation patterns, and polynya events are discussed in this chapter. 
Following this background discussion. Chapter III delves into the physics of the coupled 
ice-mixed layer system, examining both previous scientific studies and the current work. 
Requirements for the accurate modeling of this physics are also discussed. Chapter IV 
reviews sensitivity studies, as the strengths and weaknesses of this model are discovered. 
Finally, Chapter V provides a summary of this work and details the conclusions drawn as 
well as recommendation for work that needs to be continued. 
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II. POLAR REGIONS 


This section is designed to provide the background knowledge that is important 
for the air-sea-ice model. The focus of the circulation section is on the Arctic since the 
flow patterns in the Northern Hemisphere are much more complex than the circulation 
patterns in the Southern Hemisphere. It is, however, important to note that the Weddell 
Polynya in Antarctica appeared to migrate west from year to year during its presence in 
the 1970s. Despite this fact, for the purposes of the one-dimensional vertical model, it is 
assumed that if the model can handle the horizontal advective effects of the Arctic, it is 
robust enough to handle similar effects in the Antarctic. In regards to climate, similar 
assumptions cannot be made since the atmospheric forcing can be significantly more 
extreme in the Antarctic than in the Arctic. The section concerning leads and polynyas 
focuses primarily on the Weddell Polynya in the Antarctic. The reason for this is not 
because physics does not support polynyas in the Arctic, but simply because the 
relatively large amount of information about the Weddell Polynya makes it a good 
example of and laboratory for the physics that create ice-free regimes in an ice-covered 
region. 

A, CIRCULATION OF THE ARCTIC OCEAN 

To better understand one-dimensional vertical effects in the Arctic Ocean, the 
horizontal movement of water into and through the region, as well as the bathymetry of 
the region need to be discussed (Fig. 1). For the purposes of this paper the rather 
complex circulation patterns and bottom topography of the Arctic Ocean examined by 
Rudels et al. (1994) is simplified to include two primary basins, the Eurasian Basin and 
the Canadian Basin, the Lomonosov Ridge, and the flow of water between and around 
these features. This simplification is justified since Rudels et al. (1994) demonstrate a 
strong front between the two basins at the Lomonosov Ridge. Similarly, Aagaard and 
Carmack (1994) show evidence to support the presence of a strong boundary current 
around the Lomonosov Ridge. 
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Figure 1. The Aretie Oeean with bathymetry. From Barry et al, 1993. 
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Although water can enter and leave the Arctic Ocean from several points, the 
Bering Strait will be considered the primary point of entry for Pacific water per Aagaard 
and Carmack (1994). Likewise, the Fram Strait was considered the main source for 
Atlantic and Intermediate waters according to Rudels et al. (1994). A secondary source 
of water into the Arctic is river runoff that can act to cap a water column and prevent heat 
exchange with the environment (Mellor and Hakkinen, 1994). 

Despite their proximity to one another, the two basins are different. Rudels et al. 
(1994) suggest that the circulation through the Canadian Basin is slower than that of the 
Eurasian Basin, causing the water to reside in the Canadian Basin by as much as a decade 
longer than in the Eurasian Basin. Rudels et al. (1994) show the Canadian Basin to be 
well mixed according to temperature and salinity by water type. Rudels et al. (1994) also 
depict the Eurasian Basin as being well mixed according to temperature and salinity by 
water type. However, the Atlantic and intermediate waters of the Eurasian Basin increase 
in both temperature and salinity as depth increases, whereas the same type of water in the 
Canadian Basin decreases in both temperature and salinity as depth increases (Rudels et 
al., 1994). This difference is important to the vertical heat flux and convection of each 
basin respectively. Similarly, the difference is important when water from both basin s 
come together at the Eomonosov Ridge because double diffusion may occur (Rudels et 
al., 1994). 

The circulation time is shorter in the Eurasian Basin than the Canadian Basin 
because advection is strong enough in the Eurasian Basin to pass the incoming water 
back out of the Eram Strait (Rudels et al., 1994). That does not mean that vertical 
processes do not occur in this basin. In fact, the temperature and salinity profiles detailed 
by Rudels et al. (1994) could allow entrainment of warmer underlying water into the 
colder surface layer. If this were to occur, the rising heat could melt ice. At the same 
time, the loss of heat at depth would cool the lower layers. 

The temperature and salinity structure of the Canadian Basin is of particular 
interest since Rudels et al. (1994) suggest that the structure will allow heat to be mixed 
down into the water column (Eig. 2). Mixing of heat down into the water column would 
deepen the mixed layer. 
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Figure 2. Representative temperature and salinity profiles from various loeations in the 
Northern Hemisphere high latitudes. From Coaehman and Aagaard, 1974. 
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If this mixing of heat down is greater than the entrainment up from the abyss,* or 
if the basin is not deep enough for an abyss to exist, then little heat would be added to the 
mixed layer to melt iee. If, however, iee forms beeause entrained heat is less than heat 
lost to the atmosphere, the mixed layer salinity inereases. Inereased salinity would then 
make the water sink eonvectively. This proeess ereates the warm and salty water at depth 
diseussed by Aagaard and Carmaek (1994). 

Introdueing adveetive proeesses like river runoff or narrow eurrents like the West 
Spitsbergen Current, as mentioned by Mellor and Hakkinen (1994), probably do little to 
offset the vertical processes in the Canadian Basin. River runoff would create fresh water 
at the surface that would quickly freeze and add to the ice. Currents at either the surface 
or some depth could create eddies and cause stirring in the horizontal. However, due to 
the stronger mixing of heat down in the column, these eddies would have little effect on 
the direction of vertical mixing. 

Adveetive processes like friction, wind stress, and currents undoubtedly keep the 
ice in the Canadian Basin from growing unbounded. Similarly, Lemke et al. (1990) point 
out that a deep mixed layer can resist ice growth. This situation would be very similar to 
the circumstances surrounding the presence of a polynya. However, strong mixing of 
heat down in the Canadian Basin could theoretically make it a one-dimensional ice 
generator, especially when the temperature of the water at the surface is at or near the 
freezing point. As it is, the Canadian Basin becomes a density capacitor, building up 
density over longer periods of time and then slowly but constantly discharging salty and 
relatively warm water at the Lomonosov Ridge. This water then mixes with cooler and 
fresher water released from the Eurasian Basin. When the waters from the two basins 
come in contact at the Lomonosov Ridge double diffusion may occur as noted above. 
Strong boundary currents along the Lomonosov Ridge and the east coast of Greenland 
then push and carry the water before it exits from the Fram Strait at depth (Aagaard and 
Carmack, 1994). Stable water, defined as cold salty water beneath warm fresher water, is 
the ultimate result of the double diffusion. In this way warm, salty, cold, and fresh waters 


* The term abyss is borrowed from Hakkinen (1987) and is used here to mean the depth of the water 
column beneath the mixed layer. 
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are all brought into the Aretie Oeean, split and proeessed by both the topography and 
eireulation, and then re-circulated to replenish the world’s oceans (Fig. 3). 


B, ATMOSPHERIC CONDITIONS 

Atmospheric forcing plays a role in both the northern and southern hemispheres. 
Ultimately, ice will form or melt as a result of the overall energy balance. This balance 
includes radiation fluxes, sensible heat, and latent heat (Fig. 4) in both polar regions. 
Whether open areas like leads exist, when heat is brought to the surface, is critical to the 
energy balance since Vihma (1995) states that ice cover prevents heat exchange. Okada 
and Yamanouchi (1995) similarly show that when the sea ice concentration increases in 
the Antarctic winter the surface heat flux decreases. If a large section of ice in the South 
Pole were instantaneously removed, evidence presented by Wu et al, (1996) suggests 
that the energy balance would restore the ice within one year. This not only demonstrates 
the strength of the energy balance but also suggests a seasonality to ice formation. 

Seasonal change is also an important factor in the overall energy balance. Sea ice 
concentration is greatest in the southern hemisphere in the July to August time frame 
during the Antarctic winter (Fig. 5). Similarly, the sea ice extent is greatest in the 
northern hemisphere from February to March during the Arctic winter (Fig. 6). Seasonal 
differences can be witnessed, and melting and freezing can be observed through surface 
albedo. A low albedo, measured in percentage, corresponds to melting ice and a high 
albedo classifies solid ice (Barry et al, 1993). 

There are also significant seasonal differences in ice coverage within each of the 
polar regions. Preller et al (2002) describe ice coverage in Arctic regions to be thinnest 
during the summer; from June to September, and thickest in the winter; from March to 
April. Similarly, Preller et al (2002) suggest that different seas, within the Arctic, for a 
given season will have different ice thickness and coverage, but that no sea will 
perpetually have the greatest coverage. Preller et al (2002) compare this phenomenon to 
a seesaw where for a given summer sea “A” has thicker ice than sea “B”, but the 
following summer sea “B” will have greater thickness than sea “A”. 
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Figure 3. Circulation of the Atlantic Layer and intermediate depth waters of the Arctic 

Ocean. From Rudels et a/., 1994. 
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Figure 4. Typical annual surface energy budget for the Southern Hemisphere. From Okada 

and Yamanouchi, 1995. 


10 






, , SEA ICE CONCENTRATION (60S-70S) 

K 7^ ) 

100 

60 

60 

40 

20 

0 

Figure 5. Typical annual Southern Hemisphere sea ice concentration. From Okada and 

Yamanouchi, 1995. 
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Figure 6. 



Typical annual Northern Hemisphere sea ice extent. From Walsh and Johnson, 

1979. 
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Understanding both seasonal and interannual seesaw effeets are important for the United 
States Navy when attempting to foreeast iee thiekness and extent. To do this, the Navy 
uses the Polar lee Predietion System (PIPS), (Preller et al, 2002). This model makes use 
of real time and elimatologieal atmospherie data, the Cox Oeean Model, and the Hibler 
lee model to produee a foreeast for iee approximately five days in the future (Preller et 
al, 2002) (Fig.7). One problem with PIPS doeumented by Preller et al (2002) is that it 
often over-prediets iee in the Barents Sea and under-prediets iee in the Labrador Sea. 

This over-predietion of iee may be attributed to the poor resolution of smaller adveetive 
proeesses like eurrents or river run-off (Mellor and Hakkinen, 1994). Another possibility 
for the diserepaneies in predietion is that the model does not resolve the heat fluxes in 
these regions as aeeurately as needed. Speeifieally, the Labrador Sea has heat fluxes and 
deep mixed layers more eonsistent with the Antaretie region than other Aretie seas, and 
so the vertieal proeesses of this sea generate greater amounts of iee than predieted (Guest, 
personal oommunieation). 


C. LEADS AND POLYNYAS 

Leads and polynyas are both areas of open water where iee would have been 
predieted given the elimatology of the surfaee foreing by the atmosphere. There are two 
types of polynyas enumerated by Smith et al (1990); the first is a “latent heat” polynya 
that is maintained by horizontal adveetive heating proeesses attributable to wind and 
eurrents. The seeond is a “sensible heat” polynya that is ereated from vertieal eonveetive 
heating from the oeean. In either ease the net heat flux into that area is suffieient to 
inhibit iee formation. As this researeh is foeused primarily on vertieal effeets, the 
“sensible heat” polynya is of greatest interest. Smith et al (1990) deseribe leads as open 
ehannels within the iee that are ereated by a divergenee in the iee flow. Leads are also of 
interest espeeially when eonsidering that oeean dynamies and heating from below may 
develop weaknesses in the iee. When the weaker iee is put under stress from wind or 
surrounding iee it diverges into open water ehannels. 
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Polar Ice Prediction System 2.0 



Figure 7. U.S. Navy PIPS Model. From Preller et al, 
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While lead location does have a stochastic nature, polynyas often recur in the 
same location with some periodicity. However, it is important to note, as Smith et al. 
(1990) does, that the Weddell Sea sensible heat polynya disappeared in the 1980’s after 
having appeared regularly in the late 1970’s. This indicates the strong dependence of 
open water formation on the dynamics occurring beneath the surface and suggests that 
climatological conditions, while important, are not the sole causes. Similarly, the 
seasonal nature of the heat flux, with greater heat flux out of the water in winter, is no 
assurance of ice formation; nor does summer warming of the water column guarantee 
complete ice melting. 

In regards to the Weddell Sea in Antarctica, Martinson et al. (1981) assert that the 
Weddell Polynya cannot be created by wind forcing or atmospheric effects alone since 
the low-pressure system required to create it would not survive long enough in this 
environment. Instead, they postulate that some type of preconditioning occurs to raise 
the pycnocline from below the mixed layer. Martinson et al. (1981) offer several 
possibilities for the source of the preconditioning including both topography and mean 
flow spin-up. The preconditioning, whatever its source, coupled with salt rejection at the 
surface during freezing causes a layer of cold fresh water above a raised portion of warm 
salty water from below. When the density of the cold fresh water surpasses that of the 
warm salty water the column can overturn leaving warm salty water at the surface. In 
this way freezing of ice initiates a negative feedback where the end result is ice melting 
from the heat of the newly raised water just below the ice at the surface of the water 
column. Because of density driven overturning Martinson et al. (1981) contend that the 
temperature and salinity profiles of both the pycnocline beneath the mixed layer and the 
mixed layer itself are vital for determining the formation of the Weddell Polynya. It is of 
interest to note, as well, that although no conclusions are made in regards to the effects of 
topography, Martinson et al. (1981) do speculate that the polynya is subject to the 
circulation of the Weddell gyre and appears to cycle from the east to the west with the 
mean flow. Furthermore, they suggest evidence that the polynya moves to an extreme 
westward location before it is recreated at its extreme eastern-most point (Martinson et 
al, 1981). 
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III. THEORY 


A. PREVIOUS STUDIES 

Recognizing the importance of atmospheric conditions and seasonal variations is 
necessary for predicting ice freezing and melting. Similarly, as in the case of the 
Weddell Polynya and mentioned by Martinson et al. (1981), the circulation patterns, both 
surface and subsurface, and/or bathymetry may aid or precondition melting or freezing 
effects. However, since low or high pressure systems and their accompanying winds 
come and go at varying time scales, and bathymetry remains relatively constant within a 
region, the physics of the water column is the most significant component of sustained 
surface freezing and melting. More specifically, the coupled thermodynamics and 
dynamics that together mix water vertically upwards and downwards contribute the most 
to ice freezing or melting. Lemke et al. (1990) explain the importance of this coupling by 
suggesting that a situation may develop where ice melts due to thermodynamic effects 
creating an area of no ice, but then freezes again, almost immediately, due to water 
movement at the surface. This argument highlights the need for the source of the heat, in 
addition to the heat it provides, to be brought to the surface for a sustained no-ice area to 
exist in the polar climate. 

Hakkinen (1987) explains the thermodynamics and dynamics necessary to bring 
heat to the surface and sustain it using a two-layer model. The two layers suggested by 
Hakkinen (1987) are a cold fresh surface layer and a warm salty bottom layer that is at 
rest. As described by Kraus and Turner (1967) the heat and salinity fluxes between these 
two layers are dependent upon the difference between the temperature and salinity 
profiles of each layer respectively. Hakkinen (1987) investigates the vertical transfer of 
heat on the mesoscale. Specifically, she focuses on entrainment of warm salty water 
from the lower layer up into the surface layer at the ice edge in the marginal ice zone 
(Hakkinen, 1987). Her one-dimensional model suggests that entrainment has an inverse 
dependence on density; as the depth of the mixed layer increases entrainment decreases, 
thereby ensuring that the mixed layer does not deepen infinitely (Hakkinen, 1987). 
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Dynamically, Hakkinen (1987) relies on strong winds blown aeross the surfaee to 
ereate drag on the water eolumn. This drag, whieh is three times as great when there is 
iee at the surfaee as when there is no ice present, creates strong Ekman transport 
(Hakkinen, 1987). Not surprisingly, the greater the drag the greater the Ekman transport, 
and so Ekman transport is greatest under the iee (Hakkinen, 1987). In order to quantify 
the drag Hakkinen (1987) defines a fraetional ealeulation of iee eompaetness, A. 
Hakkinen’s (1987) preeonditioning, as termed by Martinson et al. (1981) for the priming 
of upeoming interaetion between the two layers, eomes from the stronger Ekman 
transport down beneath the iee that leads to inereased upwelling at the iee edge. Through 
this meehanism upwelling assists entrainment in keeping warmer and more saline water 
at the surface and keeps the area free of ice (Hakkinen, 1987). 

Eemke et al. (1990) similarly use a one-dimensional two-layer model eonsistent 
with Kraus and Turner (1967) to explain sustained heat at the surfaee. Eike Hakkinen, 
Eemke et al. (1990) is eoncerned with the momentum balanee between the atmosphere, 
iee, and water column. Eemke et al. (1990) eombines the equations of motion, iee 
eompaetness and eoverage, the nonlinear viseous-plate constitutive law from Hibler 
(1979), and atmospherie foreing to explain vertical movement of heat. Eemke et al. 
(1990) also find that the greatest vertieal heat flux in the Weddell Sea is at the iee edge. 

Unlike Hakkinen (1987), however, Eemke et al. (1990) foeus less on atmospherie 
forees like wind and more on the entrainment veloeity as the meehanism for heating at 
the surface. In other words, the overall oeeanic heat flux is dependent on the entrainment 
veloeity (Eemke et al., 1990). This means that whether ice melts or freezes it is directly 
dependent on the entrainment veloeity. Aeeording to Eemke et al. (1990) positive 
entrainment brings heat to the surfaee and melts the iee, entrainment less than zero 
extends the mixed layer to the depth of the water eolumn, and entrainment equal to zero, 
setting the mixed layer depth equal to the Mon in-Obukhov length, oecurs when the two 
layers are thermodynamieally independent of one another. Eemke et al. (1990) formulate 
the entrainment values to be dependent upon the equation of state, and thereby effeetively 
emphasize the role of the temperature and salinity expansion eoeffieients. Eurthermore, 
the inelusion of the equation of state underscores the importance of ehanges in both 
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temperature and salinity between the abyss and the mixed layer for vertieal heat exehange 
(Lemke et al, 1990). 


B, THERMOBARICITY 

Thermobarieity is an important proeess to inelude when attempting to model or 
understand the eomplete physies of the mixed layer. It is defined as the eombined 
dependenee of seawater thermal expansion on temperature and pressure (MeDougall, 
1984). The work done by Lemke et al. (1990) lays the groundwork for the signifieanee 
of thermobarieity. It is important beeause it magnifies the buoyaney flux assoeiated with 
mixed layer eonveetion (Garwood, 1991). If thermobarieity is ineluded in the equation of 
state the thermal expansion eoeffieient, a , is no longer a eonstant, but is instead a 
funetion of both temperature and depth. This means that entrainment is not only a 
funetion of the differenees in temperature and salinity between the abyss and mixed layer, 
but is also a funetion of the depth and temperature of the abyss-mixed layer boundary. 

The physieal result of the inelusion of thermobarieity is eonsistent with Lemke et al. 
(1990); when entrained heat from below is greater than the heat leaving the water 
eolumn, the mixed layer warms and any existing iee begins to melt. Similarly, if 
entrained heat is less than the heat leaving the eolumn the mixed layer eools and freezing 
oeeurs at the surfaee. However, the magnitude of the entrainment eould be very different 
with the inelusion of thermobarieity, under certain conditions, resulting in either freezing 
or melting when Lemke et al. (1990) predict the opposite. 

Garwood et al. (1994) shows the necessity of understanding thermobarieity and 
the importance of being able to model it. He enumerates four different scenarios where 
thermobaric effects have implications in polar modeling. 

1) Mixed-layer turbulent kinetic energy (TKE) may be increased several-fold by 
the thermobaric amplification of the buoyancy flux (Garwood, 1991). This phenomenon 
is not significant at temperate and low latitudes, but it becomes large for polar sea mixed 
layers. It cannot be neglected for realistic air-sea ice interactions. Large eddy simulation 
(LES) verifies its importance in increasing entrainment. However, commonly used 
parameterizations of ocean mixing, including second-order closure of Mellor and 
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Yamada (1982) and k-profile parameterization (KPP) (Large et al., 1994), do not yet 
inelude this physies. 

2) Buoyancy is not conserved in coastal subduction. Although buoyancy is well 
conserved in vertical advection of water masses at lower latitudes, it is not a conservative 
property in the polar seas. This is especially important for vertical circulation that is 
enhanced on shelf slopes. Here it becomes a severe limitation for layered models that 
presume a constant potential density for each layer (Garwood and Isakari, 1993; Jiang 
and Garwood, 1995, 1996, 1998). 

3) "Thermobaric stability" is an increase in the hydrostatic stability of seawater 
having positive spiciness and nearly density-compensating T-S structure. It causes spicy 
warmer and saltier parcels to resist vertical mixing because of a newly recognized 
thermobaric component to the buoyancy frequency. This added component depends 
upon the thermobaric depth . This phenomenon enables internal gravity waves 

without vertical buoyancy gradient (Garwood, 2002), and it helps explain lateral intrusion 
of filaments of spicy water across T-S frontal zones. Furthermore, the theory predicts 
reduced stability for cold and fresh filaments with negative spiciness. This process is 
potentially important in enhancing shoreward transport but reducing the flux of 
suspended and dissolved material from the shelf into the Arctic interior. 

4) A critical mixed-layer depth for maintenance of open water and sensible 

heat polynyas (Smith et al, 1990) under freezing atmospheric temperatures is predicted 
because of thermobaricity (Garwood, 2001). Whenever the mixed layer depth exceeds a 
critical depth, the entrainment heat flux will exceed the surface heat loss to the 

atmosphere , resulting in ice melting regardless of the rate of heat loss to the 

atmosphere from the ice-mixed layer system. This is a positive feedback mechanism for 
warming the surface layer, limiting ice formation, and cooling the ocean interior. It 
undoubtedly plays a role in establishing the extent and geographic evolution of the 
marginal ice zone, and in the maintenance of polynyas and their entrainment of nutrient- 
rich water into the surface layer. 
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C. FUNDAMENTAL EQUATIONS 

Since Kraus and Turner (1967) first modeled the mixed layer, oeeanographers 
have worked to understand the role of turbulenee in mixed layer dynamies. Their foeus 
was limited to the temperate oeeans and the role of air-sea interaetions on turbulenee. 
Independently, since Hibler (1979), polar oeeanographers have proposed and advanced 
models for iee mechanics, with cursory attention to mixed layer turbulenee. Recognizing 
the need to eonsider mixed layer dynamies under the ice, Lemke et al. (1990) and 
Hakkinen (1987) ineluded linear formulations for entrainment into their coupled iee- 
oeean numerieal models. However, Garwood (1991) showed that the equation of state’s 
nonlinear properties needed to be ineluded when eonsidering turbulent mixing in the eold 
water of the polar seas. 

The physieal proeesses of the eoupled system depend heavily on the vertieal 
entrainment of warmer, saltier, and denser water from below. Normally in the lower 
latitudes the denser water from below is both eolder and saltier than the water above it, 
but in the high latitudes where salinity is the dominate determinate of density, it is the 
warm and saline water that makes up the high density water from below. Entrainment, or 
the rate of deepening of the mixed layer, is important for two primary reasons. First, 
entrainment is needed to ealculate both the temperature and salinity fluxes that make up 
the buoyaney flux. Seeond, when ice is present at the surfaee, entrained water may heat 
the ice layer directly without heating the mixed layer significantly, and so ice thiekness is 
directly dependent on the entrainment. 

Prediction of polar sea mixed layer dynamies and thermodynamies is dependent 
upon the same physies that governs all water pareel motion: the Reynolds equations for 
the mean momentum, which are derived from the Navier-Stokes equations, 

U,+UU^+Wy+WU^=- — + fv - fyW-{u'u''^ -{u'v'^ -{u'w'^ +v{uxx+Uyy+Uzz^{\) 
U^+UU^+Wy+WU^=- — +fv - fyW-{u'u'^ -{u'v'^ +v{uxx+Uyy+Uzz^{T} 

w^+uw^+^y+ww^ =- — +fyil-g -(vV) -(wV) +v{wxx + Wyy+Wzz^{l>) 
the continuity equation. 
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(4) 


u,+Vy+w^=0 

the heat and salt equations where the Reynolds deeomposition is also applied with 
T = T + T',S = S + S', averaged, and put into flux form as, 

!^ = T, + uT, + W^+wT^ - (^), - (^'), + ’^r(T„+f„ + t_) (5) 

= S, +ii5, + W, + ivJ = -(SV)^-(SV')^ -(SV')^ + *•,(S„ + + SJ (6) 

and the equation of state, 

p = pM-a(T-Tj + Ji(S-Sj] (7) 


Predietion of the turbulent fluxes, and entrainment in partieular, is based upon the 
Turbulent Kinetie Energy (TKE) equation, whieh is also derived from the Navier-Stokes 
equations. 
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This equation ealeulates the energy available to do work ereated from shear produetion, 

, and the vertieal fluxes of temperature and salinity that make up the 


,du 

dz 


U'w' - h V V- 


dz 


vertieal buoyaney flux, b'w' = agT'w'- /3gS'w'. The four terms to the right of the 
equals sign in the TKE equation are, from left to right, as follows; the shear produetion of 
horizontal turbulent kinetie energy, the produetion or damping of vertieal turbulent 
kinetie energy by buoyaney forees, redistribution, and dissipation. These terms apply to 
the mixed layer whether or not iee is present at the surfaee. The shear produetion term is 
always non-negative, and horizontal turbulent kinetie energy is a direet result of the wind 
stress at the surfaee and produees “foreed eonveetion.” Vertieal turbulent kinetie energy 
is produeed whenever the buoyaney flux is positive due to negative salinity flux and/or 
positive temperature flux. Positive buoyaney flux results in “free eonveetion.” 

Buoyaney flux will always be negative at the bottom of the mixed layer, where the 
turbulent kinetie energy is used to entrain the underlying heavier water into the mixed 


22 



layer. Redistribution moves turbulent kinetie energy vertically within the mixed layer. 
Dissipation converts turbulent kinetic energy to heat, and is always a loss term. 


To solve for the entrainment it is necessary to solve for the vertical integral of the 
TKE equation. If equation (8) is vertically integrated over the mixed layer, the 
entrainment velocity may be calculated in terms of the stratification (AT and A^), mixed 
layer depth, h, and the forcing velocity scales (w, and w,), which are defined below. 


w„ 
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a^ghAT 
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The resulting equation (9) shows that all of the available energy to do work, or stirring in 
this case, is provided from either forced convection from wind or ice forcing, or from the 
free convection from within the water column. Forced convection. 


CjM, 
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dz dz 
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originates in the atmosphere and quantifies the turbulence created from wind blowing 
across the surface or across ice at the surface. Free convection, or positive c^wi ', is the 
vertical integral of the buoyancy flux over the mixed layer, 

= f b'w'dz (11) 

J -h 

As can be seen in the equation for entrainment (9) the equation of state (7) has a 
significant influence on the buoyancy flux and magnitude of free convection (11). It 
should be mentioned that the term “free convection” only applies when wj’ is positive. If 
is negative the integrated buoyancy flux dampens the turbulence. Appropriately, this 
situation is referred to as “buoyant damping.” The value of can be positive under the 
ice if freezing is occurring at the surface. If melting is occurring wj’ is negative 


regardless of the value. In the ice-free system is positive when the net heat flux 
is positive (out of the system). When the net heat flux is negative, or into the ice-free 
system (summer conditions), is negative as well. 
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The net effect of dissipation is absorbed by the coefficients, c^ and Cj. For this 
work the coefficients Cj and Cj are considered tunable constants. Kraus and Turner 
(1967) assumed q and Cj to be constant. First Gill and Turner (1976), and later 
Garwood (1977) showed that both q and Cj vary with mixed layer stability. However, 
Cj and Cj may be considered constant for slowly changing surface forcing conditions, 

which are assumed in this study. The parameter Ha is the "thermobaric depth" 
(Garwood, 1991) for which increased pressure doubles the thermal expansion relative to 
the surface value, , and will be discussed further in the following paragraphs. 

The buoyancy flux is a function of both the potential temperature flux and the 
salinity flux. The vertical integral of the buoyancy flux, 

0 _ 0 _ _ 

J b'w\z)dz = J {b'Wg' + b'Wg'^dz 

—h —h 


shows this dual dependence. Breaking buoyancy flux into these two components and 
solving demonstrates how each of the fluxes is dependent on vertical entrainment. The 
two components are. 


{\)b'w,' = ag 
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The vertical fluxes are related to entrainment by. 


(1) O'w'-h = -tS.9w^ 

( 2 ) '¥w'-h=-^Sw^ 


where IS.0 and are the jumps in potential temperature and salinity at the base of the 
mixed layer. The potential temperature (^ ) and the in situ temperature (T) are assumed 
to be equal in this study, with negligible error. Taking the vertical integral of each 
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component yields , the vertieal integral of the buoyancy flux, and its dependence on 
the vertieal entrainment veloeity: 


= - PoSh 


S-S, 


(^ogh* 1 + — 


P.gh*^S 

2 


The temperature and salinity fluxes follow directly from the temperature and salinity 
profiles as well as the net heat flux into or out of the water column (Fig. 8). Temperature 
and salinity profdes in the high latitudes are quite unique, especially with the presenee of 
ice (Fig 9). Assuming ice eoverage, the water just below the ice will be warmer than the 
iee; if it were not warmer it too would be frozen. Consequently, the water just below the 
ice is warmer than the freezing point and usually continues to get warmer with depth. 
Similarly, since ice rejects salt when it freezes, the ice will have a much lower salinity 
than the water beneath it. Salinity usually continues to increase with depth. Reealling 
that salinity dominates over temperature in the high latitudes, these profiles create a 
stable density structure. 


When there is ice coverage all of the entrainment works direetly to heat the ice 
layer and does not heat the mixed layer signifieantly. When this occurs the ice begins to 
melt and so the vertical entrainment velocity affects ice thickness (hi) as well; 

= = 03) 

dt Pi^f Pi^f 

where mf is the melting-freezing term, quantifies the latent heat of fusion, and pc^ 
is the water density times the heat capacity. When the entrainment heat flux exceeds the 
upward flux of heat and radiation at the ice surface (Qq), then the ice will melt making 
dk 

< 0. Conversely, when the entrainment is weaker than the surface heat flux iee will 


freeze making —- > 0 . 

dt 
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Figure 8. Representative temperature and salinity flux profiles for the high latitude regions. 
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Atmosphere 



Figure 9. Representative salinity and temperature profiles for the high latitude regions. 
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In this way vertical entrainment simultaneously affects both free convection and the ice 
thickness. Therefore, to solve for entrainment requires solving three simultaneous 

dh- 

equations, one for free convection, wj ', one for change in ice thickness,- - , and one 

dt 

for vertical entrainment velocity itself, . 

Accuracy of the calculation of buoyancy flux is improved by considering 
thermobaric enhancement of entrainment. This occurs when the thermal expansion 
coefficient, a, is no longer treated as a constant but is allowed to vary with both 
temperature and depth. Consequently, the equations for both vertical entrainment 
velocity and free convection contain the term/Za, the "thermobaric depth", as well as the 
surface term for the thermal expansion coefficient, . The inclusion of thermobaric 

depth can lead to a magnified entrainment rate, increased buoyancy, and more efficient 
stirring of the mixed layer through free convection. It is through this process that heating 
of the ice can occur more quickly and with greater intensity. Prolonged enhanced 
entrainment in a region can lead to a polynya, where the ice has been melted away 
completely and the entrainment works to warm the mixed layer to greater depth. 

Under certain situations free convection can be turned off, = 0 , or even be 

made negative. As discussed previously, when this term is negative it is no longer called 
free convection but is instead called buoyant damping. This situation occurs when the 
entrainment shuts down or is very weak. Entrainment, however, can never be made 
negative since un-stirring, or directly reversing completed work cannot be achieved. If 
free convection were to shut down, forced convection could dominate and wind stirring 
of the mixed layer could result. When this happens the mixed layer could shallow to the 
Obukhov Length (L), and a step would be noticed in the density profile consistent with 
summer shallowing conditions. It is important to note that for the model created in 
association with this work, circumstances leading to a negative sum of the free and forced 
convection will not be processed accurately. 


D, MODEL 
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A one-dimensional model eonsisting of a system of ordinary differential equations 
was solved as an initial-value problem numerieally to examine the interaetion between 
the mixed layer and iee at the surfaee. The model was designed this way so that 
partieular attention eould be paid to the importanee of thermobarieity as well as the 
strength of various foreing terms and initial eonditions. Additionally, this model was 
built with the idea that it eould be used direetly or as the preeursor for a model that will 
be able to fit with eurrent iee predietion models. Speeilieally, the model should 
ultimately work in eonjunetion with the Hibler lee Model so that the United States Navy 
ean use it operationally in systems like PIPS (Fig .7). The end goal for the future model 
is simple; foreeasts for iee formation, thiekness, and extent must all improve over today’s 
model output. For example, the iee thiekness in the Greenland Sea, whieh is eurrently 
often poorly predieted, may be more aeeurately represented. The model must be robust 
enough to generate iee-free areas like polynyas or iee eovered regions given any initial 
iee eondition and representative atmospherie eonditions. In other words, given realistie 
foreing and initial eonditions, the model should prediet a polynya in an area that is iee 
eovered or freezing in a polynya. Matlab’s C++ was ehosen as the programming 
language for this model beeause of its matrix, differential equation, and graphieal 
eapabilities, as well as its ability to perform large ealeulations with minimal user effort. 

1, Model Equations Summarized 

a. Coupled Mixed Layer-Ice System Equations 

The system of eight equations is derived by vertieally integrating the 
budgets for enthalpy (heat)(5), momentum (2)(3), salinity (6) and TKE (8) between the 
iee surfaee and the base of the turbulent mixed layer. Drag (momentum flux) and heat 
flux between the iee and mixed layer are ealeulated with bulk aerodynamie sealing laws, 
with the temperature of the iee-water eontaet assumed to be at the freezing point, Tfp. 
Entrainment fluxes of heat, salinity and momentum are derived from integrating the 
respeetive eonservation laws (2, 3, 5, and 6) aeross the entrainment zone diseontinuity 
(Garwood, 1977). 

Mixed Layer Depth; 
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dh I 3 
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Mixed-Layer Temperature; 
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Mixed-Layer Salinity: 
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Mixed-Layer Currents; 
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Ice velocity; 
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Momentum and Heat Flux Between Ice and Mixed Layer: 


T.-PCD[ur^)4{ui-u) +(v/-v) ( 22 ) 
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2. Numerical Methods 

Matlab is a programming language that ean perform large calculations and 
provide graphical representation of data through the manipulation of matrices and 
execution of user produced or existing functions. Its ability to quickly and accurately 
perform complicated and or lengthy computations make it an attractive tool for ocean 
modeling. 

a. Use of Stiff Solver 

ODE23s is a Matlab function that solves ordinary differential equations 
using a stiff solver. The use of ODE23s allowed for the tracking and graphing of many 
different oceanographic parameters as the program stepped forward in time. Specifically, 
it kept track of and depicted how both the mixed layer depth and ice thickness increased 
or decreased given the initial and forcing conditions. The program ODE23s was chosen 
over other Matlab differential equation solvers because of its ability to solve stiff 
differential equations. A way of solving this type of equation became necessary when 
approaching or crossing into the ice-free regime, where intermediate solutions change 
rapidly, but the final solution is found on a longer time scale (Shampine, 1994). Not 
using a stiff solver would mean allowing for the time step to decrease so that every 
intermediate solution could be calculated. This ability was not necessary since the final 
solution was the solution of interest. To numerically solve these types of equations 
Matlab employs the Runge-Kutta formula or, in the case of ODE23s, a modified 
Rosenbrock formula that is more efficient (Math Works, 2002). Additionally, because of 
the computational power and time benefits of this program, it was verified that treating 
the water column as three slabs was a robust method for the desired modeling. This was 
proven graphically by showing that the total of the velocities of the ice layer slab and 
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mixed layer slab, individually, were equal in value to the veloeity of a combined ice- 
mixed layer slab. 
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IV. SENSITIVITY STUDIES 


The following sensitivity studies simulated Weddell Sea winter conditions. 
Akitomo (1999) first applied Garwood et al. (1994) thermobaric instabilities to the 
Weddell Sea, but without the ice-mixed layer feedback model used here. This was done 
because there is greater interaction between temperature and salinity in this region and at 
this time. For each simulation one of the forcing conditions, net heat flux (2^) or wind 

speed, was held constant while the other was allowed to vary. Similarly, the initial 
conditions were varied for the simulations. Two different depths of the mixed layer were 
chosen to examine the effects of the ratio between mixed layer depth and thermobaric 
depth. These depths were 100 m for a moderately deep mixed layer and 1000 m for a 
deep mixed layer. The effects of changing the ratio of thermal stratification to salinity 
stratification, a^gAT to /^ggAS , on ice melting and freezing at the surface were 
examined. 

It should be noted that the program allows for the ratio of CC^gAT to P^gAS to 
change with changing temperature and salinity values. If these values were instead 
prescribed, accuracy would be lost since it would effectively not allow for any storage of 
heat or salt in the mixed layer when ice is present at the surface. Furthermore, any warm 
salty water entrained would not change the mixed layer temperature and salinity values. 
Equations from the previous chapter demonstrate how important the temperature and 
salinity fluxes within the mixed layer are to the free convection and total turbulent kinetic 
energy of the system. Not allowing for these fluxes is the equivalent of neglecting 
physical processes essential to the coupled ice-mixed layer system. Instead, by allowing 
for this change in the ratio, the vertical integral of the buoyancy flux can be accurately 
calculated. 

In all simulations the initial ice thickness is 10 cm; the initial temperature is the 
temperature of the freezing point at the surface, -1.66 Celsius; the initial salinity value is 
34.52 parts per thousand; and the thermobaric depth is 93.6 m. Since both 100 m and 
1000 m mixed layers are deeper than the thermobaric depth, thermobaricity is expected to 
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enhance the entrainment of warmer and saltier water from below. Each case is run for 
ten days, or 8.64x10^ seconds. 

A, MIXED LAYER DEPTH EQUAL TO 100 METERS 

For the first simulation at h=100 m, the ratio of cx^gAT to ^^gAS equaled 0.45 for 

the surface temperature and salinity values. The total heat flux was a constant 0.001 
cal/cm s upward and out of the coupled system. The wind speed was allowed to vary 
from 400 cm/s to 1600 cm/s in 400 cm/s increments. The results show that for each wind 
speed positive entrainment deepens the mixed layer (Fig. 10). The weakest wind, 400 
cm/s, produced the least amount of mixed layer deepening. The resultant entrainment 
was not strong enough to melt ice at the surface and instead fostered ice growth. 

The mixed layer salinity increases for the 400 cm/s wind as freezing occurs and 
salt is rejected from the ice. The corresponding mixed layer temperature warms until it 
comes into equilibrium. Once the temperature equilibrium is reached, the rate of salinity 
increase, ice growth, and mixed layer deepening all become constant. 

The higher wind values of 800 cm/s, 1200 cm/s, and 1600 cm/s each produced 
strong enough entrainment to melt the ice and deepen the mixed layer, with the 1600 
cm/s wind having the strongest melting rate. Similarly, the 1600 cm/s case had the 
largest decrease in salinity associated with the melting of the ice and injection of fresh 
water into the mixed layer. Correspondingly, the mixed layer temperature warmed the 
most significantly immediately before eventually approaching equilibrium. It should be 
noted that the 1600 cm/s plots do not extend to the full ten days like the other wind 
speeds. This is due to the strong entrainment that occurs and the ill-defined matrix that 
results when trying to continue to solve for the entrainment and terms 
simultaneously. 


34 



c 


ppt 


-1.56 


-1.58 


- 1.6 


-1.62 


-1.64 


- 1.66 


Mixed Layer Temperature vs. Time x lO”* ''S-Time 



400cm/s 

800cm/s 

1.15 



— 1200cm/: 




1600cm/E 

1.1 





cm 



1.05 



1 



0 2 4 6 8 10 


0 2 4 6 8 10 


t/sec X 10 = 

Mixed Layer Saiinity vs. Time 


t/sec 

ice Thickness vs. Time 


x 10^ 


34.56 

34.54 

34.52 

34.5 

34.48 



10 






5 

- 


cm 



0 

■ 





-5 

' 

, . , , 

-10 

. 


0 2 4 6 

t/sec 


8 10 
X 10^ 


0 2 4 6 8 10 

t/sec X 10® 


Figure 10. Simulation #1: Sensitivity study for varying winds with h =100 m, Qo =0.0010 
eal/cm^s, and ratio of Gf^gAT to /3QgAS equal to 0.45. Wind speed ehanges from 400 
ern/s to 1600 cm/s in 400 cm/s intervals. 
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Winds of 1200 cm/s and 1600 cm/s are capable of melting away the iee 
eompletely; however, onee the iee is melted away the model shows negative iee thiekness 
at approximately 2.3 days (2.0x10^ see) and 6.9 days (6.0x10^see) respeetively. This is 
an unrealistie depietion of the true physieal proeess; there is no sueh thing as negative iee. 
Similarly, an unrealistie deerease in salinity values is aehieved for these wind speeds at 
the same points in time. The salinity values are unrealistie beeause onee the iee has 
melted away there is no other souree of fresh water to eontinue to deerease the salinity of 
the mixed layer. Future versions of the model need to be able to melt away the iee, 
remain iee-free without predieting negative iee, and then be able to make iee, beginning 
with no iee, when eonditions warrant iee growth. For this situation the iee should melt 
away due to strong entrainment and then the system should remain iee free due to a 
sustained high value of entrainment. The high entrainment rate should eontinue 
deepening even without iee at the surfaee, although the true deepening rate is probably 
different than the one represented here with negative iee thickness. 

Examining the plot of mixed layer temperature with time reveals that there is 
more heat storage in the mixed layer beneath freezing eonditions (wind speed equals 400 
em/s) at the surfaee than there is with melting eonditions (800-1600 em/s wind speeds) at 
the surfaee. That is, the mixed layer equilibrium temperature is warmer beneath iee 
growing eonditions than it is beneath iee melting eonditions. This makes sense when the 
iee at the surfaee is thought of as a eap that does not allow the entrained heat to eseape 
out of the system. Sinee it is not strong enough to melt the iee and eannot eseape to the 
atmosphere, this heat is stored in the mixed layer. Also shown through this sensitivity 
study is that for a given net heat flux, wind ean eause entrainment to be strong enough to 
melt iee even in freezing atmospherie eonditions. Furthermore, strong entrainment ean 
eause and sustain polynya-like eonditions. 

For the seeond simulation at h=100 m, the ratio of (X^gAT to P^gAS equaled 
0.036 for the surfaee temperature and salinity values. The total heat flux was a eonstant 
0.001 eal/em s upward and out of the eoupled system. The wind speed was allowed to 
vary from 400 em/s to 1600 em/s in 400 em/s inerements. The results of this wind 
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sensitivity study show how important the ratio of cx^gAT to J^^gAS is to entrainment and 
iee formation or melting. 

Like the first simulation, the mixed layer depth inereases for all wind speeds due 
to positive entrainment (Fig.l 1). The 1600 cm/s wind caused stronger entrainment than 
the 400 cm/s wind as well. However, for this case the entrainment is never strong enough 
to melt the ice. Entrainment strength is still dependent on the wind speed, with the higher 
wind speeds creating greater entrainment. This can be seen in the freezing rate where the 
lower wind speeds with weaker entrainment have higher freezing rates. In this case the 
400 cm/s wind freezes at a rate of approximately 1.1 cm/day (11 cm in ten days) as 
compared to the 1600 cm/s wind that freezes at a rate of approximately 0.4 cm/day (4 cm 
in ten days). 

Examining the ice thickness time series plot shows that the 1600 cm/s wind slows 
the freezing rate. Given a few more days with this constant wind speed it is conceivable 
that the ice formation would halt and that melting might even begin. Salinity predictably 
increases as freezing occurs at all wind speeds. Here, however, the salinity increases the 
most for the 1600 cm/s wind as compared to lesser winds. This is in contrast to the result 
from the first simulation where the greatest freezing rate had the greatest increase in 
mixed layer salinity. Comparing the 400 cm/s wind to the 1600 cm/s wind and the 
resulting ice formation of both it is seen that the 400 cm/s wind produces a weaker 
entrainment, greater ice thickness, and a less saline mixed layer as compared to the 1600 
cm/s wind that produces greater entrainment, thinner ice, and a more saline mixed layer. 
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Figure 11. Simulation #2; Sensitivity study for varying winds with h =100 m, Qo =0.0010 
cal/em^s, and ratio of a^gAT to ^^gAS equal to 0.036. Wind speed changes from 400 

cm/s to 1600 cm/s in 400 cm/s intervals. 
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The 400 cm/s wind is the slowest to bring the system into thermal equilibrium, 
but appears to achieve a higher equilibrium temperature than some of the greater wind 
speeds. This can be seen by looking at the mixed layer temperature plot at day ten where 
the 400 cm/s wind has already surpassed the mixed layer temperature value for the 800 
cm/s wind and appears close to exceeding the mixed layer temperature of the 1200 cm/s 
wind. Like in the first simulation, where ice is present and freezing is occurring, the 
mixed layer has a higher equilibrium temperature leading to more storage of heat from 
the effect of ice capping the system. These results indicate that strong entrainment with 
ice presence will bring the system into thermal equilibrium faster, but that entrainment is 
not the only factor for establishing what that point of thermal equilibrium will be. 

For the third simulation at h=100 m, the ratio of a^gAT to /^^gAS again equaled 

0.45 for the surface temperature and salinity values. The wind speed was a constant 
value of 700 cm/s. The heat flux was allowed to vary from 0.0005 cakcm s to 0.0025 
cal/cm s in 0.0005 cal/cm s increments. The results show that for all values, 

sufficient entrainment was produced to deepen the mixed layer (Fig. 12). Similar to the 
second simulation, the stronger the heat flux the greater entrainment it caused and the 
deeper the mixed layer grew. A low Qq value equal to 0.0005 cal/cm s causes deepening 

to approximately 103 m, while a value equal to 0.0025 cal/cm s causes deepening to 
115 m, for a difference of 12 meters of deepening between the two fluxes. 

For all cases of the ice actually grows at first but then begins to melt. The 
time that the transition from freezing to melting occurs differs for each value; 
however, the greater the Qq value the longer it takes for this transition to occur. For 
example, for a heat flux, gg equal to 0.0005 cal/cm s it takes approximately 1.2 days 
(1.0x10 sec) for melting to begin. A heat flux, , equal to 0.0025 cal/cm s will take 4.0 
days (3.5x10^ sec) for this change from freezing to melting to occur. 
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Figure 12. Simulation #3; Sensitivity study for heat fluxes with h =100 m, wind speed =700 
cm/s, and ratio of CCf^gAT to /^QgAS equal to 0.45. Heat flux changes from 0.0005 
cal/cm^s to 0.0025 cal/cm^s in 0.0005 cal/cm^s intervals. 
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Mixed layer salinity inereases when ice is growing and decreases when ice is 
melting. This is consistent with the previous two simulations and the physics of salt 
rejection during ice freezing. When equal to 0.0005 cal/cm s begins melting ice, the 
salinity predictably starts to decrease, since the melting puts fresh water into the mixed 
layer. However, at larger values of when melting occurs the salinity does not 

decrease as it does with the equal to 0.0005 cal/cm s; instead, the rate of salinity 
increase becomes constant. The explanation for this is that entrainment is lifting more 
salinity into the mixed layer from below than it is taking out of the mixed layer through 
melting and fresh water addition to the mixed layer. 

High heat flux values, like 0.0025 cal/cm s, that lead to relatively strong freezing 
had the highest mixed layer equilibrium temperature, and also approached equilibrium 
faster than lower heat flux values. These results can be observed by comparing equal 
to 0.0005 cal/cm^s with equal to 0.0020 cal/cm^s. For 0.0005 cal/cm^s the mixed 
layer temperature came to equilibrium at approximately -1.64°C after 3.5 days 
(3.0x10 sec), while the 0.0020 cal/cm s heat flux came to equilibrium at approximately - 
1.597C after 8.7 days (7.5x10^ sec). For these measurements the rate at which the mixed 
layer temperature comes to equilibrium for the lower heat flux is .0057C/day and the rate 
at which the mixed layer temperature comes to equilibrium for the higher valued heat 
flux is .0069C/day. These results suggest that there may be a relationship between ice 
thickness and heat storage. This relationship would come into practice when strong 
entrainment does little melting but mixed layer depth continues to increase (deepen). 
Here, where the ice is not melted and heat cannot escape to the atmosphere, the heat 
becomes stored in the mixed layer and causes warming of the equilibrium temperature. 

Comparing this simulation to the first simulation, where the heat flux was kept 
constant, produces some interesting results. First, the lowest value produces the 

coolest equilibrium temperature, whereas the lowest wind speed appeared to have the 
warmest equilibrium temperature. This indicates that more heat can be stored in the 
mixed layer of the first simulation as compared to the mixed layer of the third simulation 
for minimum forcing conditions. Second, a 400 cm/s wind with constant produced 
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2 

freezing of ice, but the lowest value, 0.0005 cal/cm s, with constant wind gives 

melting. Both of these observations can be explained by concluding that the constant 
values used in simulation are weak enough to allow other factors to influence the 
entrainment and thus change the characteristics of the coupled system. The usefulness of 
these constant values for heat flux and wind speed is further substantiated because in 
simulation 1 where gg is constant and winds are weak, freezing is forecast. Similarly, in 

simulation 3 where wind speed is constant and is weak, melting is forecast. 

For the fourth simulation at h=100 m, the ratio of to /^^gAS equaled 0.036 

for the surface temperature and salinity values. The wind speed was a constant value of 
700 cm/s. The heat flux was allowed to vary from 0.0005 cal/cm s to 0.0025 cal/cm s in 
0.0005 cal/cm s increments. The results demonstrate that positive entrainment causes 
mixed layer deepening for all heat flux values (Fig. 13). The ice thickness growth rate 
increases as Qq increases from 0.0005 cal/cm s to 0.0025 cal/cm s with no transition to 
melting. Salinity predictably increases for all heat flux values, with the strongest rate of 
salinity increase corresponding to equal to 0.0025 cal/cm s. Like simulation 3, the 

greatest change in ice thickness is associated with the warmest equilibrium temperature 
and the largest amount of heat storage within the mixed layer. 
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Figure 13. Simulation #4: Sensitivity study for heat fluxes with h =100 m, wind speed =700 
cm/s, and ratio of a^gAT to /3f^gAS equal to 0.036. Fleat flux changes from 0.0005 
cal/cm^s to 0.0025 cal/cm^s in 0.0005 cal/cm^s intervals. 
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When comparing this simulation with the previous one the effects of the different 
ratios of cx^gAT to ^ogAS can be examined. It appears as though the smaller ratio, 
0.036, increases the significance of while simultaneously weakening the role of 
entrainment. A weakening of entrainment would explain why the mixed layer does not 
deepen as much in simulation 4 as it does in simulation 3. The strongest value for Qq in 
the third simulation deepens the mixed layer to approximately 115 m. The strongest 
value for in the fourth simulation only deepens the mixed layer to 103.5 m. Similarly, 

the strongest value of Qq in the third simulation produces only approximately 11.5 cm of 
ice, while the strongest value of in the fourth simulation produces approximately 37 

cm of ice. The fact that the mixed layer temperature remains cooler for each heat flux 
value in simulation 4 as compared to the same heat flux values in simulation 3 lends 
credence to the possibility that entrainment has weakened, since less heat is stored in the 
mixed layer of simulation 4. Further evidence of weakened entrainment is found in the 
mixed layer salinity versus time plot where for simulation 4 the salinity values for each 
heat flux value reach higher values than they do in the third simulation. This corresponds 
to the net effect of increased salt rejection and no addition of fresh water from melting 
caused by entrained heat. 

B, MIXED LAYER DEPTH EQUAL TO 1000 METERS 

For the first simulation at h=1000 m, the ratio of a^gAT to P^gAS equaled 0.45 
for the surface temperature and salinity values. The total heat flux was a constant 0.001 
cal/cm s up and out of the coupled system. The wind speed was allowed to vary from 
400 cm/s to 1600 cm/s in 400 cm/s increments. The results show that positive 
entrainment deepens the mixed layer for all of the wind speeds (Fig. 14). The 400 cm/s 
wind causes sustained deepening by creating relatively weak entrainment that is not 
strong enough to melt the ice. The three other larger values for wind cause freezing 
initially but then shift to melting conditions when entrainment overpowers the net heat 
flux. This initial freezing period lasts approximately 2.3 days (2.0x10^ sec). After this 
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period, the mixed layer deepens twiee as fast in the 400 cm/s case as it does in the 800 
cm/s case. 

This rate can be determined by picking values off of the mixed layer depth versus 
time plot for both the 400 cm/s line and 800 cm/s line corresponding to times of 
2.0x10^ seconds and 4.0x10^ seconds (4.6 days). For the 400 cm/s wind at time equal to 
2.0x10^ seconds the mixed layer depth equals 1013 m, and at time equal to 
4.0x10^ seconds the mixed layer depth equals 1021 m. For the 800 cm/s wind at time 
equal to 2.0x10^ seconds the mixed layer depth equals 1010 m, and at time equal to 
4.0x10^ seconds the mixed layer depth equals 1014 m. Solving for the change in depth 
per day for each case yields a change of 3.47 m/day for the 400 cm/s wind and a change 
of 1.74 m/day for the 800 cm/s wind. 

The ice thickness versus time plot changes according to the strength of 
entrainment where greater entrainment causes melting to occur sooner. Predictably, 
salinity increases with increasing freezing and decreases for strong melting, as is the case 
with the 1600 cm/s wind. The 400 cm/s wind has the highest equilibrium temperature of 
the four wind speeds. Conversely, the 1600 cm/s wind has the lowest equilibrium 
temperature but also achieves thermal equilibrium more quickly than the lesser wind 
speeds. This result is due to the fact that the high winds produce strong entrainment that 
is enhanced by thermobaricity to melt the ice without significantly warming the mixed 
layer. Like thermobaricity, the ratio of mixed layer depth to thermobaric depth likely 
plays a role in the sustained deepening of the mixed layer for the 400 cm/s wind. The 
ration of the mixed layer to the thermobaric depth was calculated using depth of the 
mixed layer, equal to 1000 m, divided by thermobaric depth, which is equal to 93.6 m for 
all of the simulations. 
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Figure 14. Simulation #5: Sensitivity study for varying winds with h =1000 m, Qo =0.0010 
cal/em^s, and ratio of Gf^gAT to /3QgAS equal to 0.45. Wind speed changes from 400 
cm/s to 1600 cm/s in 400 cm/s intervals. 
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For the second simulation at h=1000 m, the ratio of (X^gAT to P^gAS equaled 

0.036 cal/cm s for the surface temperature and salinity values. The total heat flux was a 
constant 0.001 cal/cm^s up and out of the coupled system. The wind speed was allowed 
to vary from 400 cm/s to 1600 cm/s in 400 cm/s increments. The results show that 
positive entrainment causes deepening for all cases, where the stronger the wind was the 
greater the mixed layer deepened (Fig. 15). For this simulation entrainment never 
became strong enough for melting to occur as evidenced by the constant rates of 
deepening for all wind speeds. The freezing rate was strongest for the lightest wind but 
all of the wind speeds had very similar freezing rates. Not surprisingly increase in 
salinity for all wind speeds was observed. 

Looking closely at the temperature scale on the y-axis of the mixed layer 
temperature versus time plot shows that the mixed layer temperature began near 
equilibrium and remained relatively constant for this time period. Since there is very 
little heat storage in the mixed layer with ice present at the surface it seems possible that 
entrainment is either very weak or has not yet become significant as compared to the net 
heat flux. 

As compared to the previous simulation this simulation’s mixed layer does not 
deepen as much for any of the wind speeds. Similarly, the mixed layer remains cooler 
and very little heat storage takes place for all the wind speeds in this simulation as 
compared to the fifth simulation. On average the mixed layer salinity increases more for 
the sixth simulation than it does for the fifth simulation. However, for the 400 cm/s case 
the salinity only reaches 34.525 ppt as compared to the slightly more saline value of 
34.526 ppt from simulation 5. 
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Figure 15. Simulation #6; Sensitivity study for varying winds with h =1000 m, Qo =0.0010 
cal/em^s, and ratio of a^g/ST to ^^gAS equal to 0.036. Wind speed ehanges from 400 

ern/s to 1600 cm/s in 400 cm/s intervals. 
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For the third simulation at h=1000 m, the ratio of (XffgAT to J^^gAS equaled 0.45 for the 

surfaee temperature and salinity values. The wind speed was a eonstant value of 700 
em/s. The heat flux was allowed to vary from 0.0005 eahem^s to 0.0025 eal/em^s in 
0.0005 eal/em s inerements. The results show that positive entrainment eauses deepening 
of the mixed layer for all values of Qq (Fig. 16). There is a greater entrainment for the 
higher values of heat flux demonstrated by the inereased deepening of the mixed layer for 
the higher values of heat flux. For example, for equal to 0.0025 eal/em s the mixed 

layer depth inereases to approximately 1055 m, whereas the mixed layer depth for 
equal to 0.0005 eal/em s only inereases to approximately 1010 m. For all values of 
entrainment eauses a shift from freezing to melting eonditions. The heat flux equal to 
0.0005 eal/em^s transitions first from freezing to melting at day 4.6 (4.0x10^ see). Eaeh 
subsequent heat flux transitions later than its predeeessor until gg equal to 0.0025 

eal/em^s transitions at day 5.2 (4.5x10^ see). This transition oeeurs more quiekly for 
lower heat flux values beeause even though the entrainment produeed by this low heat 
flux value is less than the entrainment produeed by higher heat flux values, it is more 
easily overeome by the entrainment. 

Salinity follows predietably from the mixed layer depth plot. Again, even though 
melting oeeurs the salinity eontinues to inerease due to entrainment of higher salinity 
water from below. Signifieant heat storage in the mixed layer ean be seen where the 
mixed layer equilibrium temperature is warmest for the highest value. As before, the 

highest equilibrium temperature appears linked to the highest iee thiekness so that 
entrainment below iee brings heat up into the mixed layer and holds it there while the 
entrained heat works to melt the iee. 
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Figure 16. Simulation #7: Sensitivity study for heat fluxes with h =1000 m, wind speed =700 
cm/s, and ratio of (Xf^gAT to /^QgAS equal to 0.45. Heat flux changes from 0.0005 
cal/cm^s to 0.0025 cal/cm^s in 0.0005 cal/cm^s intervals. 
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Comparing this simulation with simulation 5 so that the effects of keeping gg 

constant versus keeping wind speed constant produce interesting results. For this 
simulation the mixed layer depth can increase by as much as 25 m over the maximum 
mixed layer depth of simulation 5. However, in simulation 5 the maximum depth was 
produced by the 400 cm/s wind, whereas in this simulation the deepest mixed layer was 
produced by the 0.0025 caFcm s heat flux. This demonstrates that the constant values for 
gg and wind speed employed are weak enough to allow other factors to influence the 
entrainment and thus change the characteristics of the coupled system. 

Comparing this simulation with the third simulation shows the effects of differing 
initial mixed layers. The effects include; increased change of mixed layer depth from 
approximately 15 m for the greatest heat flux in simulation 3, to 55 m for the strongest 
heat flux in simulation 7; increased maximum salinity value for the greatest heat flux in 
simulation 3 (34.547 ppt), as compared to the maximum salinity value for the greatest 
heat flux in simulation 7 (34.530 ppt); and mixed layer equilibrium temperature while 
relatively constant for all of the heat flux values between the two simulations does 
increase faster in simulation 3 than in simulation 7, where equilibrium is reached in 
simulation 3 after approximately 3.5 days (3.0x10^ sec) as compared to 8.7 days 
(7.5x10^ sec) for simulation 7. 

For the fourth simulation at h=1000 m, the ratio of (X^gAT to P^gAS equaled 

0.036 for the surface temperature and salinity values. The wind speed was a constant 
value of 700 cm/s. The heat flux was allowed to vary from 0.0005 caFcm s to 0.0025 
cal/cm^s in 0.0005 cal/cm^s increments. The results show that positive entrainment 
creates mixed layer deepening for all heat flux values (Fig. 17). Entrainment is weak 
since the mixed layer does not increase much with time (maximum of 3.5 m), freezing 
occurs for all heat flux cases with the maximum ice thickness equal to 39 cm for 

equal to 0.0025 cal/cm s, and very little heat can be stored in the mixed layer since the 
mixed layer temperature varies little even in the most extreme heat flux case (0.003C 
over 10 days). 
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Comparing this simulation with simulation 4 (h=100 m) to examine the effeets of 
initial mixed layer depth produees the following: 3.5 meters of depth inerease for both 
simulations where equals to 0.0025 eal/em s; simulation 8 produees approximately 2 

em more iee for all heat flux eases than simulation 4; the mixed layer temperature warms 
faster and to larger equilibrium values in simulation 4 than it does in simulation 8; and 
the mixed layer salinity inereases faster and to higher values in simulation 4 than it does 
in simulation 8. These results ean be verified by looking at the slope of the respeetive 
lines or by taking sample points off of the lines and eomparing the ehange in a quality 
over time numerieally as done on several oeeasions above. 

Comparing this simulation with simulation 6 to examine the effeets of gg versus 

wind foreing on the deep mixed layer indieates that the effeet of strong heat flux values 
out of the system is more signifieant to freezing than low wind speeds. This ean be 
quantified by ealeulating the iee growth from simulation 6 for low wind (400 ern/s) to be 
equal to 12 em. Similarly, ealeulating the iee growth from simulation 8 for high heat flux 
(0.0025 eal/em s) yields 30 em of iee. From these results it ean be eoneluded that for 
these values of heat flux and wind speed, large heat flux tends to produee stronger 
freezing eonditions than low wind foreing. 
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Mixed Layer Temperature \s. Time 
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Figure 17. Simulation #8; Sensitivity study for heat fluxes with h =1000 m, wind speed =700 
cm/s, and ratio of a^gAT to /^^gAS equal to 0.036. Heat flux changes from 0.0005 
cal/cm^s to 0.0025 cal/cm^s in 0.0005 cal/cm^s intervals. 
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Comparing this simulation with simulation 7 to examine the effeets of the ratio of 
a^gAT to /3(^gAS shows that the higher ratio, like that of simulation 7, produees melting 

of the iee and greater deepening of the mixed layer (up to 55 m gained as eompared to a 
maximum of 32 m gained in simulation 8). Simulation 7 also ean have greater heat 
storage within the mixed layer as a result of the inereased warming of the mixed layer. 
The mixed layer salinity plots, while reaehing similar values for the time length of the 
simulations, ehange at different rates. The rate for simulation 8 is nearly linear for all 
heat flux values and is equal to .0013ppt/day for equal to 0.0025 eal/em^s. The rate 
for simulation 7 is not linear sinee the entrainment for this simulation is strong enough to 
eause a transition between freezing and melting eonditions. This rate is equal to 
O.OOlppt/day for equal to 0.0025 eal/em s. While this differenee in rates appears 
small, examination of the plots shows that for simulation 7 the rate appears to tend 
toward equilibrium and may, after a long enough period, beeome negative or deerease the 
salinity of the mixed layer. Meanwhile, the rates for the various heat fluxes in simulation 
8 eonstantly inerease the salinity of the mixed layer. Ultimately, eomparing the results of 
simulation 8 with simulation 7 suggest that the low ratio of (X^gAT to ^^gAS from 
simulation 8 eannot support as strong an entrainment rate as the higher ratio of CC^gAT to 
Pag AS from simulation 7. 
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V. CONCLUSIONS 


A. SUMMARY 

The coupled ice-mixed layer system is a delicate balance of physical processes 
that are dependent on both initial conditions and forcing conditions. Changes in the heat 
and or buoyancy fluxes, wind speed, initial ice thickness, and initial mixed layer depth 
can all initiate feedback processes that can lead to either ice growth or melting. This 
balance of physical processes is so tenuous that a moderate change in any quantity or 
combination of quantities rarely produces an intuitive result. Based on this balance it is 
easy to understand why different basins in the Arctic at constant latitude can have 
different amounts of ice thickness at the surface. Similarly, polynya formation and 
maintenance is better understood in light of these intricate relationships. Ultimately, ice 
formation or melting cannot be forecast from surface forcing conditions alone, as has 
been attempted in the past. Thorough knowledge of the site-specific temperature and 
salinity profiles and fluxes are necessary in conjunction with the surface forcing 
conditions to predict entrainment and ice thickness. 

Thermobaricity is a real and significant factor for quantifying the strength of 
entrainment and the turbulent kinetic energy of the coupled system. In comparing 
simulations between different initial mixed layer depths the thermobaricity seems 
particularly important or sensitive to wind forcing as compared with changing heat 
fluxes. In other words, the greatest changes in mixed layer depths were found for the 
varying wind cases with an initial mixed layer depth of 1000 m. For these cases the 
change in mixed layer depth could be as large as 60 m. The comparable case for initial 
mixed layer depth of 100 m could only achieve a maximum of 15 m. (To see these 
maximum possible depths compare Figures 12 and 16.) This scenario is especially 
important for sustaining polynya events and indicates that perhaps surface wind forcing 
does more to maintain polynya conditions than net surface heat flux. Furthermore, this 
mechanism for rapid and strong deepening of the mixed layer is important because it 
likely helps explain deep-water regeneration on a global scale. 
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Another significant result of the simulations was an observed relation between 
greater ice thickness and higher heat storage in the mixed layer. In particular as 
described above every simulation shows that for the greatest freezing rate it takes the 
mixed layer temperature longer to come into equilibrium and thus stores and holds on to 
more heat as it comes into equilibrium. Simulation 2 and simulation 6 (in Figures 11 and 
15, respectively) appear to be exceptions to the rule, however. For their given conditions 
it appears that the simulation did not run long enough to allow the mixed layer 
temperature to come into equilibrium. This can be seen where the plots are still rising up 
and starting to bend toward equilibrium. This result is significant because it indicates 
that thick ice at the surface acts as a cap on the coupled system. This cap does not stop 
entrainment from melting the ice at the surface; instead, melting begins before the system 
comes to thermal equilibrium indicating that when there is ice at the surface the entrained 
heat works to melt the ice before it works to warm the mixed layer. 

Perhaps the most significant result of the simulations was to verify the complexity 
of the interactions that occur in the coupled system and which are represented in the 
equations developed in the theory section. For example, as predicted by the equations, no 
one factor was determined to be the most important to quantify for determining 
entrainment strength. Instead, because of the feedbacks inherent in the system, a 
condition that led to strong entrainment in one case may not lead to strong entrainment in 
another case, depending on the other initial and forcing conditions. Because of these 
feedbacks, it is critical that modeling of coupled ice-mixed layer system include the 
realistic thermodynamics described above. Similarly, this work has demonstrated that the 
PIPS model would become significantly more useful for ice prediction if the necessary 
physics were included. 


B, RECOMMENDATIONS FOR FURTHER STUDY 

There are many possible topics of further study for high latitude mixed layer 
dynamics and thermodynamics. Preconditioning for polynyas by advective processes at 
depth like warm subsurface eddies, as mentioned by both Hakkinen (1987) and Lemke et 
al, (1990), should be explored as a possible mechanism to aid entrainment in polynya 
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initiation and maintenance. From an atmospheric perspective, further research of the 
effects of precipitation in ice-free and ice-covered regions is important for greater 
resolution of heat flux data and the total energy balance. Net surface heat flux into the 
system needs to be thoroughly explored just as net heat flux out of the coupled system 
was examined in this work. Similar sensitivity studies of net heat flux both into and out 
of the coupled system need to be conducted for an initial ice thickness value equal to 
zero. For this study it is critical that the model be able to transition freely and correctly 
between the ice-free system and the ice-covered system. 
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